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Abstract 

Massive Compact Objects (Machos) are currently being discovered at substan- 
tially higher rates than would be expected from standard models of known stellar 
populations. To determine whether they are due to non-standard distributions of 
known populations or to a heretofore unknown ('dark') population, one must ac- 
quire more information about the individual events. Space-based parallaxes are po- 
tentially the best tool for extracting additional information, yielding the "reduced" 
speed V = v/{l — z), the "reduced" Einstein radius fg = [4:GMDqq,z/ (1 — z)]^^"^ /c, 
and the transverse direction of motion $. Here M is the Macho mass, v is its trans- 
verse speed, and z is its distance relative to the source distance, Dqq. For example 
bulge and disk Machos can be cleanly separated by measuring v. To leading order, 
parallax measurements by a single satellite result in a four-fold degeneracy: two 
possible values of v and two possible signs for the component of motion perpen- 
dicular to the projected satellite-Earth vector. It had been believed that a second 
satellite would be required to break this degeneracy. I show that the velocity dif- 
ference between the satellite and the Earth allows one to partially or totally break 
the degeneracy using a single satellite. Eor most Macho events it is possible to 
measure v and fg. Eor some it is also possible to measure $. The proposed Space 
Infrared Telescope could measure ~ 100 parallaxes per year by applying ~ 400 hr 
of telescope time. 
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1. Introduction 



Several dozen candidate lensing events have now been detected toward the 
galactic bulge by the OGLE (Udalski et al. 1994 and references therein) and 
MACHO (Alcock et al. 1994; C. Stubbs, private communication 1994) collabo- 
rations. Estimates of the observed optical depth to lensing in this direction of 
T ~ 3-4x10"^ far exceed the sum of the values expected for a standard disk, 
~ 5 X 10~^ (Paczyhski 1991; Griest et al. 1991) and a standard axisymmetric 
bulge, T;, ~ 5 X 10~^ (Kiraga & Paczyhski 1994). Whether the additional observed 
events are due to additional material in the disk, an elongated bulge, a highly 
flattened halo, or to some other structure, they are potentially very exciting. Ex- 
trapolating from current detection rates, I estimate that the MACHO collaboration 
could see ~ 100 events next bulge season if they monitor 35 fields. 

The MACHO (Alcock et al. 1993; K. Griest, private communication 1994) and 
EROS (Aubourg et al. 1993) collaborations have seen a total of 5 candidate events 
toward the Large Magellanic Cloud (LMC). Although far fewer than expected 
from a standard spherical halo, this is also far more than would be expected from 
a combination of the known stars in the disk, thick disk, and spheroid (Gould, 
Miralda-Escude, k Bahcall 1994; Bahcall et al. 1994) or from the LMC itself 
(Gould 1994c). Thus the source of these events is also a puzzle. 

In the course of detecting a Massive Compact Object (Macho) by means of its 
time-varying magnification of a source star, one generally obtains only one piece 
of information about the Macho itself, the time scale lo~^ . This is a combination 
of three of the four physical parameters that one would like to know: 

vc 

where M is the Macho mass, v is its transverse speed relative to the observer-source 
line of sight, and z is the ratio of its distance -Dql to that of the source -Dqs? 

(1.2) 
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The source distance -Dqs is assumed to be known from spectroscopic parallax. The 
fourth potentially measurable parameter is $, the direction of transverse motion. 
The direction does not enter the time scale lo~^ . The distribution of observed time 
scales is generally the only information available about the Machos. Even when this 
is measured over several lines of sight, the implications of such measurements are 
highly degenerate, especially because we are certainly seeing two distinct Macho 
distributions and very probably more than two. Hence, in order to understand the 
nature of the events it is necessary to extract more information about them. 

Space-based acquisition of Macho parallaxes is potentially the most powerful 
method of obtaining additional information (Gould 1992; Gould 1994b, hereafter 
Paper I). First, Macho parallaxes yield two additional pieces of information about 
each event. When combined with the time scale, the three parameters are the 
"reduced" Einstein radius fg, the "reduced" transverse speed C, 



^2 _ ^ GDqs Mz 

1 — z' 1 — z' 



fg = — -, v = - -, (1.3) 



and $, the transverse direction relative to the satellite-Earth separation vector. Eor 
example, the reduced speed can discriminate well between disk and bulge Machos 
seen toward the bulge, or between Galactic and LMC objects seen toward the 
LMC (Han & Gould, in preparation). Once these populations are distinguished, 
the reduced transverse speed and direction give direct kinematic information about 
each population, and the reduced Einstein ring gives information about the mass 
spectrum. A second important feature is that once the appropriate satellite is 
launched, space-based parallaxes can be obtained for a large fraction of all events. 
By contrast, special circumstances are required to obtain ground-based parallaxes 
(Gould et al. 1994) or proper motions (Gould 1994a, Nemiroff & Wickramasinghe 
1994; Maoz & Gould 1994). 

Parallaxes are determined by measuring the difference in the times of maximum 
magnification At, and the difference in the impact parameters A/3, as determined 
from the Macho light curves seen from the Earth and a satellite. See Eigures 1 and 
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2 of Paper I. The light curves are described by the magnification 

= /f + xit) = [uHt - hf + (1.4) 

where to is the time of maximum magnification and fi is the impact parameter in 
units of the Einstein radius, = [4:GMDoszi^~z)V^'^ /c. The vector displacement 
in the Einstein ring 

Ax = (cuAt, A/3), (1.5) 

is related to the two-dimensional separation vector r between the satellite and the 
Earth by, 

Ax = — . (1.6) 

re 

Unfortunately, while coAt = io{t'Q — to) is well defined by comparing the value of t'^ 
measured from the satellite light curve with to measured from Earth, the impact- 
parameter difference A/3 has a four fold degeneracy: A/3 can have two distinct 
magnitudes, 

A/3± = |/3'±/3|, (1.7) 

corresponding to impact parameters on the opposite or same side of the source. 
These different values for A/3 lead to corresponding different values for Ax± = 
[(cuAt)^ -|- (A/3-1-)^]"^/^ and so for v = ijjr/Ax±. In addition, A/3 can be of either 
sign. The two signs correspond to the two directions perpendicular to the projected 
satellite-Earth separation r, i.e., $ = it tan~"^(A/3/c<;At). In Paper I, I discussed 
the possibility of breaking this four-fold degeneracy by making observations from a 
second satellite. However, since the launching of even one satellite into solar orbit 
is a formidable undertaking, the need for a second satellite would be a significant 
drawback. 
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Here I show that the degeneracy in A/3 can be broken by measuring the differ- 
ence in time scaies 

Au = u'-u, (f.8) 
between the sateiiite and the Earth. 

2. Breaking the Degeneracy 

Let the three-dimensionai separation vector between the sateiiite and the Earth 
be denoted R, and let their relative velocity be denoted U. Let s be the direction 
of the source. Then the projected separation r and projected velocity u are given 
by 

r = R-(R-s)s, u = U-(U-s)s. (2.f) 

I define the parallel and perpendicular components of u by 

r • u (r X u) • s , , 

u|| = , = ^ (2.2) 

I assume that the projected velocity difference is small compared to the reduced 
speed u <^ V. The difference in inverse time scales Aco is then given by 



(2.3) 



Aio |v + u|-?; ?5||M|| + ?5_LM_L 

id V 

Since V = cufg, equation (f .6) (or eq. 2.8 from Paper I) implies, 

ru^ At ruiAfi rco 

^11 = T7 — ^-L = 71 — 775 ^ = 7 — • (2-4) 

" {Axf {Axf Ax ^ ^ 

Thus, equation (2.3) becomes 

u _i A [3 = r Aid - u\\uj At. (2.5) 

Since u and r are known from the satellite's orbit and Aco and coAt are known 
from the comparison of light curves, A/3 is unambiguously determined by equation 
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(2.5). The only question is: is A/3 well enough determined to break the four-fold 
degeneracy which arises in the previous analysis? 

For definiteness I will assume that (5^, the error in estimating Alo/uj from the 
two light curves is given by 

6^ = 0.01, (2.6) 

and that the errors in estimating cuAt, A/3_, and A/3_|_ are respectively 

St = 0.76^, S-=0.5S^, S+ = 8S^. (2.7) 

In § 6, I will discuss the reasons for the relative error sizes given by equation (2.7) 
and the prospects for achieving the accuracy given by equation (2.6). 

If information about Aco is ignored, then the surface over the Ax = {to At, A/3) 
plane will have four minima. The minima at (cuAt, ibA/3_) will have 1 a error el- 
lipses that are 0.007 X 0.005 and the minima at (cuAt, ibA/3_|_) will have 1 a error 
ellipses that are 0.007 X 0.08. To break the degeneracy, the information encoded 
in Aid must distinguish between these minima at the ~ 3 cr level. 



3. Satellite With An Earth-Like Orbit 

To address the problem of degeneracy breaking concretely, I assume that the 
satellite is in an Earth-like orbit, but displaced from the Earth by an angle 6. Then 
R = 2sin(6'/2)i?®, U = 2 sin(6'/2)u® , and R • U = 0. Here = SOkms-^ is the 
Earth's orbital speed and i?® = 1 AU. Let a be the angle between s and the south 
ecliptic pole, and let t/j be the phase of the orbit such that at = 0, the projected 
displacement vector r is most closely aligned with s. I then find 

r = R{1- sin^ a cos^ V)^^^ u = 0®i?(l - sin^ a sin^ V)^^^ 

R o • 2 sin 2V' (3-1) 



Uj_ 



-flmR — cos a, Mil = flffiR — sin a- 



r ' w 2 ' 

where 0® = t^®/i?® = 27ryr~"^. 
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4. Observations Toward the LMC 



Since the LMC is very nearly at the south ecliptic pole, a ~ 0, equation (3.1) 
implies that r = i?, u = [/, and 

M_L = M|| = 0. (4.1) 

Hence, equation (2.5) becomes 

The smallest separation in A/3 between the four solutions will generally be 0{Ax). 
In any event, this is the smallest separation that needs to be resolved since if two 
solutions are much closer than this they represent essentially the same reduced ve- 
locity. Combining equations (2.4), (2.6) and (4.2), I find the condition for complete 
degeneracy breaking: u±_ ^ 36i_jV or 

Vm R 



Ra 



^ 0.03 (vector solution). (4-3) 



Note that for Machos in the Galactic halo, C/u® ~ 10, while for LMC Machos 
5/^0 ^ 50. From equation (4.3), it might appear that the best strategy is to 
place the satellite near 6 ~ 180° so that R ~ 2i?0. In fact, one must be careful 
that R ^ fe since otherwise the event as seen from the satellite will be out of the 
Einstein ring. Equation (4.3) shows that the degeneracy can be broken for Galactic 
Machos provided that R ^ 0.3 i?® but cannot be broken for LMC Machos unless 
R > Rq) or the accuracy is improved beyond equation (2.6). It is much easier, 
however, to distinguish between the ibA/3_ solutions from the ibA/3_|_ solutions 
since these are separated by (9(1) rather than 0{Ax). Erom equations (2.6) and 
(4.2) I find the less restrictive condition: 0® ^ 3(5^cu or 



id 



^ 1.7 days (scalar solution). (4-4) 



When equation (4.4) holds, v and fg can be measured, but $ cannot be measured 
unless equation (4.3) also holds. 
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5. Observations Toward the Bulge 

Observations toward the bulge cover many lines of sight. For definiteness, I 
adopt a line of sight with an angle a = 84° from the south ecliptic pole, similar to 
Baade's Window. That is, cos a = 0.10. 

The fact that the bulge lies near the ecliptic introduces new complications in 
the parallax measurement which are not present for the LMC. The size of the 
displacement in the Einstein ring Ax = (cuAt, A/3) between the Earth and the 
satellite is a strong function of the time of year. Since Ax = r/fe^ the measurement 
is most effective if the projected displacement is as large as possible without being 
larger than the reduced Einstein ring. Erom equation (3.1) it is clear that r varies 
by a factor sec a ~ 10 over the course of a year. Thus, no matter what physical 
separation R is chosen and no matter what is the distribution of fg, the projected 
separation will be non-optimal during some parts of the year. Nevertheless, the 
difficulties presented by the changing length of r are not in and of themselves very 
serious. Typical events arising from disk Machos have fg ~ 2 AU (Af/0.1 Af©)"^/^ 
while those in the bulge have fg ~ 5AU(M/O.lM0)i/2. During ~ 80% of the year 
the projected separation is in the range 0.33 < r/R < 1. Hence, if i? ~ i?® then 
over most of the year and for most Machos, 0.04 ^ Ax ^ 1, so that according 
to equation (2.7), it should be possible to measure Ax to within ^ 20% (up to a 
four-fold degeneracy). 

The fact that the bulge lies near the ecliptic does pose significant problems for 
breaking the degeneracy. The error in estimating A/3 in the degeneracy-breaking 
equation (4.2) has two terms, [ruj /uj_)6i_j and [u^^/uj_)6t. As I discuss below, the 
second term can effectively be ignored. Thus the condition for complete degeneracy 
breaking is u±_ ^ 36i_jV or 

V R 

— (1 -|- tan^ a sin^ r/j)""^/^ ^ 0.03 (vector solution). (5-1) 

V R(Q 

Note that this expression is smaller than the ecliptic-pole formula (4.3) by a factor 
(l-|-tan^ a sin^ . This factor remains in the range 0.1-0.2 for about 2/3 of the 
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year. For disk Machos, typical reduced speeds are v ~ 300kms~ , implying that 
for most of the year few disk lensing events can be completely resolved. Since bulge 
Machos have substantially higher C, it will generally not be possible to completely 
resolve these events. 

However, it will generally be possible to distinguish between the A/3_|_ and 
A/3_ solutions and thereby measure v and fg for the majority of Machos presently 
being detected. Following the same argument that led to equation (4.4), I find 
0® ^ 3Scoio or 

^ 1.7(1 + tan^ a sin^ r/j)"^/^ days (scalar solution). (5-2) 

I now justify ignoring the error arising from the second term in equation (2.5). 
From equations (2.7) and (3.1), one finds that the ratio of the errors of the second 
relative to the first term is 0.35 (O® /cu)(l — sin^ a cos^ V')""^ sin(2?/'). This is greater 
than unity only ii lo~^ > 17 days. Fven in this case, the second term is significant 
for only a short fraction of the event. Hence, it has no major effect. 

6. Error Estimates and Systematic Effects 

I derived the estimates of the error ratios (2.7) by assuming that a "typical" 
lensing event (/3 ~ /3' ~ 0.5) is measured frequently and with uniform accuracy 
from both the ground and space over a time interval 3lo~^, beginning when the 
source enters the Finstein ring. The absolute normalization [eq. (2.6)] was set by 
assuming that the photometry has errors of cr = 0.01 mag and is carried out with 
frequency / = lOco. For different parameters, the errors scale oc af~^^^. 

A very important implicit assumption of the calculation is that the satellite 
and ground-based observations are carried out using very similar filters. The im- 
portance of this assumption can be seen as follows. Under the observational as- 
sumptions listed above, it is typically possible to measure the impact parameter /3 
of an individual light curve only to an accuracy of ±0.12. Similarly it is possible 
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to measure u only to a fractional accuracy of ±18%. At first sight these high 
errors make the estimates for the uncertainties in A/3_ and Acu given in equations 
(2.6) and (2.7) appear ludicrous. However, the reason for the high errors in fi and 
Lo is that one must allow for light from an additional unresolved star, either the 
lensing star, a companion to the lens or the source, or a random superposition of 
an unrelated star. If such an unresolved star is present and not taken into account, 
it will make the impact parameter appear lower and the event shorter. The fit to 
the light curve is very insensitive to the fraction of light that is not being lensed, 
and errors in this estimate induce errors in fi and u. However, if the light curves 
from the ground and space are both obtained with the same filter then one knows 
a priori that the fraction of the blended source "star" due to background light is 
the same for both sets of observations. Hence the differences A/3_ and Acu can be 
determined much more accurately than the individual parameters fi and u. On 
the other hand, the sum A/3_|_ is rather poorly determined [see eq. (2.7)]. 

Let A and A' be the mean effective wavelengths of the Earth and satellite filters, 
and define AA = A — A. I then find that errors in Alo/uj as large as ~ AA/A could 
result. Thus, to avoid compromising the accuracy of the experiment, one should 
insure AA ^ 0.256^X or AA ^ 20A in / band. If this level of precision cannot be 
achieved, then the ground-based observations must be done in two bands in order 
to find the color correction term between the ground-based and space-based filters. 

There are two systematic effects that could compromise the degeneracy-breaking 
measurement of Aco: either the source or the lens could be a binary. If the source 
star suffers acceleration due to a companion with component aj_ perpendicular 
to r, the satellite and Earth light curves will differ very much as they would if 

= a_iAtz/[l — z). Eor the great majority of binary sources, this value of u±_ 
will be clearly inconsistent with being due to the known satellite-Earth velocity 
difference. On the other hand some will be consistent. Hence one must worry that 
any measured Au that is consistent with an allowed value of A/3 [see eq. (2.5)], is 
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actually produced by a companion. The inferred acceleration would be 

V cos a 1 

a± = oe — — r- 6.1 

V(s 1 — sm a cos^ tp ^; cos O 

where a® is the Earth's acceleration. In most cases this acceleration would be so 
large that it could easily be detected by spectroscopic observations (provided that 
the orbit was not very close to the plane of the sky). All sources with non-zero 
measured Aco should be monitored spectroscopically to check for possible duplicity. 

In general, the effects produced by binary lenses do not resemble those from 
satellite-Earth velocity differences so that the contamination from this effect will 
be very small. 



7. Satellite Requirements 

The optimum value of R is different for different populations of Machos. Eor 
Galactic Machos seen toward the LMC and disk Machos seen toward the bulge, 
the separation should be i? ^ i?® so that the event as seen from the satellite does 
not fall outside the Einstein ring. Eor bulge Machos and LMC Machos, R should 
be as big as possible to maximize the relatively small value of Ax. These divergent 
requirements can be partially reconciled by launching a satellite into solar orbit 
which stays in an Earth like orbit but gradually drifts away from the Earth. The 
satellite-Earth observations would then be most sensitive to nearby Machos at the 
beginning and more sensitive to distant Machos later on. 

At present, it appears that ~ 100 Machos will be discovered toward the bulge 
during each ~ 6 month bulge season. To observe each event 30 times requires 
~ 16 observations per day. If the typical star has / ~ 19, then 1% photometry 
could be obtained in 1.5 hr/day on a Im telescope or in 24 hr on a 0.25m. Similar 
accuracy could be achieved in K band assuming (/ — K)o ~ 0.5 and Ay ^ 1.5. 
A Im telescope would have the advantage of better resolution and there are many 
other imaging projects one might want to do with such an instrument. In fact. 
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the proposed Space Infrared Telescope Facility (SIRTF) could make very useful 
parallax observations with a commitment of only ~ 400hr/yr. As presently de- 
signed, SIRTF would drift 0.1 AUyr""^ behind the Farth. To optimize the satellite 
for Macho observations, this rate should be about doubled. 

Acknowledgements: I would like to thank D. DePoy for helping to clarify 
several points. 
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